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Abstract 

The main problem addressed here is to decide whether it is or not 
possible to go from a given position on a peg-solitaire board to another 
one. No non-trivial sufficient conditions are known, but tests have been 
devised to show it is not possible. We expose the way these tests work 
in a unified formalism and provide a new one which is strictly stronger 
than all previous ones. 

1 Introduction 

Peg solitaire (also called Hi-Q) is a very simple board game that appeared 
in Europe most probably at the end of the 17th century. Its prior origin is 
unknown. The first evidence is a painting by Claude-Auguste Berey of Anne 
Chabot de Rohan (1663-1709) playing it. It seems to have then become 
popular in some royal courts. The mathematical study of the game starts in 
1710 when Leibniz writes a memoir on the subject [I]. We refer the reader 
to the excellent historical account presented in Beasley's book [2]. Let us 
introduce rapidly hov^^ this game is being played. The first data is a board 
& which in first approximation may be thought of diS cl subset of Z^. The 
classical ones are the english board and the french one drawn belovi^, and we 
present a third one introduced by J.C. Wiegleb in 1779 (see [2]). 

Each square of this board can hold at most one peg, and a problem as 
we define it here is to go from a given distribution of these pegs (say /) 
to another one (say J), via a succession of legal moves that we now define. 



*keywords: Peg solitaire, Hi-Q, Pagoda function. 

tAMS classification: primary 05A99, secondary 91A46, 52B12, 90C08. 
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Figure 1: English Figure 2: French 

board board Figure 3: Wiegleb 

board 



Given three consecutive squares P, Q and i? in a row or a column (but not 
on a diagonal), of which two consecutive (say P and Q) contain a peg while 
the third one (i?) does not, a legal move consists in removing the two pegs 
in P and Q and putting one on the empty square R. We classically say that 
the closest peg in P jumps over the middle one in Q and lands in R, while 
destroying the peg in Q. As a trivial consequence, the number of pegs on 
the board decreases when the game proceeds further. For most authors, a 
problem consists in reducing the initial distribution of pegs, what we call 
thereafter the initial position, to a single peg via legal moves. They qualify 
the position as solvable if this is possible. We shall say that the problem in 
our sense is feasible if one can go from the initial position to the final one 
by using legal moves. Note that the number of such moves is known and 
equals the difference between the number of pegs in the initial position and 
the number of pegs in the final one (that is: \I\ — \ J\). 

Given a problem, we can try all possible legal moves and repeat this 
action until the required number (|/| — \J\) of moves is reached or no fur- 
ther move is possible. This process usually gets stuck because of the com- 
binatorial explosion. For instance E. Harang [3] computed that there are 
577 116 156 815 309 849 672 paths on the english board from the initial 
position consisting of the full board on which we leave the central square 
empty. Of which 40 861 647 040 079 968 lead to the final peg being on the 
central square. See also [1]. In faet, numerous setting tend to show that the 
problem is NP-complete. For this sentence to have a sense, we are to choose 
a way of extending the board to infinity, and there is no canonical fashion 
to achieve that. The case of an n x n board is studied in [5j while the k x n 
board with k fixed is shown to be linear in [6]. Of course, one may wonder 
whether the english board as a subset of a 7 x 7 board is tractable or not 
and the answer is still no, at least not without huge resources. The number 
of paths being enormous, we look for tests that will ensure us that it is not 
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possible to solve a given problem. We would welcome any test that would 
guarantee the feasibility, but none are yet known. 

The first of this test is attributed to Reiss in 1857 in [7] though Beasley 
traces it back to A. Suremain de Missery, a former officer of the French 
artillery, around 1842. We again refer the reader to [2] for more historical 
details. It is also described in Lucas book 0, which contains also more 
material and in the dedicated chapter of [9]. A seemingly more algebraic 
approach is proposed in |T0], but it turns out to be only a different setting 
for the same test. This test is very often reduced by modern authors to the 
rule-of-three test (see below). 

We shall first present these tests in a formalism that will help us clarify 
the situation; this formalism will also be adequate to present the advances 
realised on the subject in 1961/1962 at Cambridge university by a group of 
students (among which were Beasley) led by J.H. Conway. 

We shall finally present a different test, which we term quadratic, and 
which is stronger than all previous ones. It however relies on solving a larger 
integer linear program and can sometime be resource demanding. We provide 
however numerous examples that we have discovered by exploring thousands 
of problems, and this in itself shows the practicality of the approach. The 
theory of this test in its purest form is complete, but we provide in the two 
last sections several improvements of it, on which we are still working. All 
examples have been computed via an intensive use of the lp_solve library |12| . 
a GTK interface and a C-program both due to the author. 

Let us end this introduction by mentioning that Beasley also introduced 
a very geometricai tool (the in and out Theorems), but it does not fit well 
in our framework and has not been worked out for an arbitrary problem (to 
the best of my knowledge at least), even if one remains on an english board. 
We shall not discuss it here. In more recent time, there has been attempts at 
working out a model of this game via string rewriting as in [6]. This approach 
remains however fundamentally one dimensional as are string rewriting rules. 
It has had applications though in describing the complexity of the game. 

2 Main formalism of the linear board 

Given a board &, we consider the Z-module ^(©, Z) of all rational integer 
valued functions over this board, and define similarly ^((5, F2) and ^(©, Q). 
This is one of the main step of the formalization: a position in the game is 
given by a subset J C (5 (the set of squares containing a peg) , which we model 
by its characteristic function 1/. If P G (5, we note P the function that is 1 in 
P and O everywhere else. A move is thus the function f = P+Q — R and 1/ — f 



3 



should become another characteristic function; we have of course assumed 
that P, Q and R where three consecutive points in this order either in a row 
or in a column of ©. We denote the set of these moves by ^(©). In the case 
of the english board, ^((3) has cardinality 76, while (3 has cardinality 33. 

Here comes the main remark. Assume we can go from / to J by the 
succession of legal moves fi, f2, . . . , f^. Then we have 

ii-ij= J2 u- (1) 

l<i<A: 

There are three ways to exploit this writing. We can say that 

• 1/ — Ij is a rational integer linear combination of members of ^((3). 
This leads to the classical Reiss's theory, or to the lattice criterion 

of [H]. 

• 1/ — Ij is a linear combination with non-negative rational coeflficients 
of members of ^{&). This leads to the main part of Conway's group 
theory. 

• 1/ — 1 j is a linear combination with non-negative integer coeflficients 
of members of This leads to what we call the fuU linear test, or 
also the non-negative integer test. 

We introduce some notations 

v{e,z)= J2 ^-f (2) 

fe!?(s) 

and 

fGS*(6) f6®(6) 

3 Reiss theory and the rule-of-three test 

Let us first expose rapidly and in modern notations the classical material. 
Characteristic functions having values O or 1, it is tempting to look at 1/ as 
taking its values in the field with two elements F2. To avoid confusion, we 
note 1/ this characteristic function as an element of ^(©, F2). If one can go 
from the initial position / to the final one J by the succession of legal moves 
fi,f2, . . . ,fk, one still has 

i/ - ij = J2 

l<i<k 
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where fj are of course the moves seen with values in F2. U f = P + Q — R, 
then f is the function over & that takes the value 1 G F2 at all the three 
points P, Q and R, and vanishes otherwise. However, F2 is now a field and 
V{&, ¥2) is simply a vector space! Deciding whether 1/ — Ij belongs to it is 
a simple matter requiring only linear algebra. 

Let us investigate this problem further. One way to characterize V{&, ¥2) 
as a subspace of ^((5, F2) is to compute equations of it. By using the 
canonical scalar product, this reduces to computing V{@, F2)"'" which means 
the elements x ^ F2) such that 

yf = p + Q-Re^ie), x{P) + x{Q) = x{R) (4) 

since any such x verifies 

\/geV{&j2). Y,x{A)g{A)=Q. (5) 

We need a name for such elements of V{&, F2)"'", and we propose the name 
witness. Let us start to do so on the english board. Let us determine a 
function Xo- We first fix four values on a square, for instance 
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Figure 4: Starting values 
By using ([4j), we can readily extend these values: 
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Figure 5: Extension 

As it turns out, there are two ways to compute a: either by adding the 
two values on the column above its square or the two on the line containing 
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it. The result is here the same a = 1. We can use this process to compute 
the values of Xo on the full board. 

What is the dimension of V{&, ¥2) in this case? The values on the initial 
square determine the values everywhere as we have just now remarked, and 
there is thus 16 witnesses. But these values are not linearly independant and 
there are linearly generated by the four 
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We can even use this process to extend the values to Z^. This yields 



























































































































































































































































































Figure 6: Over 



Now that the reader ses the regularity of this tiling, s.he will be convinced 
that they can be extended to Z^. The way one drops the english board on it 
yields for instance this witness: 
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And once, the witnesses are determined, equations defining V{&, ¥2) are 
obtained by taking the scalar product with (a basis) of them. A classical 
problem is to determine whether it is possible to start with the french board 
filled with pegs, except for the central square that is left empty and to end 
with only one peg. This can be shown to be impossible by using the theory 
above, but we leave this pleasure to the reader. 

This theory of witnesses is essentially what is cailed Reiss's theory [7], 
though it is expressed with other words, and is present in Lucas's book [8]. 
We say "essentially" because they do not use any linear algebra and that 
their way to reach this result is by using direct move together with reversed 
ones (to undo a move). They obtain what they call characteristic positions, 
which is equivalent to the equations defining V"(©, F2). This is however what 
is presented [10]. There are still a distinction to be made: 

1. One can start from witnesses of Z^, restrict them to & and get witnesses 
for this board. This is cailed the rule-of-three. Of course, we get only a 
four independant equations that may not define V{@, F2) fully. If the 
board is thick enough, for instance when there exists a defining square 
from which all the other valnes of the witnesses can be deduced, this 
is enough. 

2. One can start from V{&, ¥2) and directly compute a basis of witnesses. 
This is required when the board is weakly connected (or even not con- 
nected!) and V{&, ¥2)'^ has dimension larger than 4. Several examples 
like that are given in [TT] . 

4 The integer linear test and the lattice crite- 
rion 

Thinking bach in terms of V{&,Z), the lattice criterion of [Ti] is to say 
that 1/ — Ij should belong to V{&, Z). How is this test connected with the 
previous one? Or, alternatively: we decided to reduce the problem modulo 2; 
Why not try to do so modulo 3? Let us first note that we may identify 
^(6, F2) with ^(6, Z)/2 ■ ^(6, Z) via 

~: ^(6,Z) ^ ^(e,Z) 

g ^ g : & ^ ¥2 (6) 
P g{P) mod 2 

During this process, V{&,Z) is of course sent on V{&,¥2)- Let us state 
formally two questions we want to answer: 
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1. Is V{&, Z) a lattice of full rank in ^(6, Z)? 

2. How to compute =F(6, Z)/V{e, Z)? 

In the sequel, we introduce a hypothesis on the geometry of the board & 
that will enables us to answer fully these questions. It will turn out that this 
will also exhibit the very tight link between the integer linear test and the 
theory of witnesses, as exposed in the previous section. 

If the two points P and R oi & are extremities of a member of ^{&), 
we say that P and R are neighbors and we note P -q- R. The reflexive and 
transitive closure of this relation is an equivalence relation, and if two points 
A and B are equivalent according to it, we note A = B. We can now state 
an important definition: 

Definition 4.1 A board & is said to be with no isolated point if for every 
point P of &, there exists a point Q = P and which is the middle point of a 
move. 

Most boards will verify this hypothesis. It means that each =-equivalence 
class contains a middle point. However the number of such classes may vary. 
For a sufficiently thick board, there will be exactly 4 classes, but there may be 
more, if the board is not connected for instance, or contains thick chambers 
very weakly connected by only one square. The reader will easily construct 
examples of boards with no isolated point but were the number of classes is 
larger than 4. The following Theorem is central in our discussion: 

Theorem 4.1 If & is with no isolated point, then 2^{&,Z) C V{&,Z). 

A final notation before sketching the proof: iff = P + Q — Re ^(6), 
we note f — —P + Q + R the reversed move (with equal middle point). 

Proof: We show that for every P G 6, we have 2P G V{&,Z). If P is 
a middle point, say of the move f, then 2P = f + f belongs to V{&,Z). 
Otherwise, there exists a chain P = Pq -©■ Pi -©■...-©■ Pn where P„ is a 
middle point. Furthermore, by definition, there exists fj G ^{&) such that 
2Pj — 2Pj+i = fi — f'i for every i = O, . . . , n — 1. Finally, we can also write 
2P„ = fn + f'n foi' some fn G Summing up all these equations, we reach 

2Po = fo-f[) + fl-f'l + --- + fn-l-f;-l + fn + f;ey(6,Z), 

which is the required conclusion since P = Pq. ooo 

This Theorem has several consequences. First of all, on such boards, 
the Q-vector spanned by the f s (that would be V{&, Q)) is the whole space: 
V{&,Z) is a sublattice of ^{&,Z) of full rank. Let us note the following 
Lemma that will be required later: 
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Lemma 4.1 If & has no is olated points, we have \&\ < \^{&)\ < 4|6| — 8. 

Proof: The lower bound comes from the faet that ^(©) generates ^{&, Q). 
For the upper bound, count horizontal and vertical moves separately. For the 
horizontal (resp. vertical) ones, count the moves according to their left-hand 
side (resp. lower) point. The lemma follows readily. ooo 

As a main consequence, we have the following Theorem. 

Theorem 4.2 Assume & to be with no isolated point and let g G ^{&, Z). 
Then 

(See ((6j) for the definition of g). 

Proof: Indeed, the direct implication is obvious, while the reversed one 
follows from Theorem 14. 11 we know that g G V{&, Z) + 2 ■ ^(6, Z) but this 
last space is nothing but V{&,Z). ooo 

This Theorem tells us that the lattice criterion is not stronger than Reiss's 
theory, when properly understood, and provided we restrict our attention to 
non-pathological boards. In faet [H] do not even give a single example when 
reduction modulo 2 does not solve the problem. Here is one: 



Figure 7: A pathological board 

The total number of pegs on the squares A, B, C and D remains constant. 
It is not difiicult to see that this example is in faet general and we have: 

Theorem 4.3 A board & is with no isolated point if and only ifV{&, Z) has 
maximal rank in ^^((3, Z). 

On boards with no isolated points, reducing the situation modulo any odd 
integer is not going to give any information; indeed Theorem 14.21 implies 
(after some work) that 

V{&, T/rnZ) = ^(6, Z/mZ) (whenever m is odd). 
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Notice finally that ^(6, Z)/y(6, Z) is simply a product of copies of Z/2Z 
in this case. It is not difficult to tackle the case with isolated points by 
generalising the reasoning used for the board drawn figure Ui and get that 
^{@,Z)/V{&,Z) is always a product of copies of Z/2Z with copies of Z. 
These results have no influence on what we develop hereafter, so we do not 
provide any formal proof. 

5 Resource counts, pagoda functions and the 
linear test in non-negative rationals 

The next main step takes place in 1961/1962 at Cambridge university when 
J.H. Conway led a group of students (among which were Beasley) that stud- 
ied this game. They came out with another and different test, also clearly 
explained in [9] and that we now describe. 

This test exploits the faet that ^ has non-negative coefiicients, i.e. the 
test consists in writing that, if we can go from I to J with legal moves, then 

Ij-ljG V+(6,Q). (7) 

As it turns out, V^{&, Q) is a cone in a vector space, and determining whether 
a point belongs to it or not is fast. We know generators of this cone (the 
elements of ^((3); they can be shown to be generator of its extreme half- 
lines), and it would be interesting to determine equations for its facets. The 
paper [T4j gives properties of these facets. In [9j as well as in |2j, so cailed 
resource counts or pagoda functions are introduced. These are functions ir 
on & such that 

Wf = P + Q-Re^{e), n{P) + n{Q)>n{R). (8) 

As a consequence, for any such function and if g belongs to l^+((5,Q), one 
has 

{7r,g) = J2A^)9iA)>0. (9) 

ags 

In particular, if one can derive J from I with legal moves, then (vr, 1/) is not 
less than (tt, Ij). Here are some examples 
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Figure 8: A re- 
source count 



Figure 9: Another 
resource count 



Figure 10: A third 
resource count 



Determining which of these corresponds to equations of facets would be 
very valuable, but their structure seems too intricate to classify them in a 
small number of regular families. For instance , a direct computation in case 
of the english board stumbles on the faet that there are an enormous quantity 
of such facets for a human eye to be able to look at them and derive some 
patterns. It is not sure that this path is blocked, though I tend to believe it 
is. 

We do not dwell any further in this part of the theory since it is extremely 
well exposed and detailed in [9], [2] and on a number of web pages. The reader 
will most probably better unterstand the strength of this theory by looking 
at section [8] of this paper. 

We should stress out here that the approach of this Cambridge group is 
commonly reduced to the use of real-valued "pagoda" functions as above. 
This is an extremely minimal understanding of their work and for instance 
does not account for the GNP balance sheet, what Beasley in [2] calls Con- 
way's balance sheet in his chapter 6; this one is however one of the main tool 
of [9]. It mixes integer valued pagoda functions together with such functions 
with valnes in F2. Beasley 's use of pagoda functions which he calls ressource 
counts (see chapter 5 of [2]) relies already on the integer character of the 
valnes taken: that is how he builds his "move map". 

The GNP diagram, or GNP balance sheet, is somewhat off our framework, 
and is in faet superseded by the next test. 

6 The linear test in non-negative integers 

The third test consists in combining both preceding ideas and write that if 
we can go from I to J with legal moves, then 



i/-ijGy+(e,z). 



(10) 
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This time, deciding that an element belongs to the integer points of a cone 
is NP-hard, but in practice, it takes only some fraction of a second on an 
english board (this was not the case in 1962!). We have of course 



V+{e,Z) c V(6,Z)p|V+(6,Q) 



fin 



and this inclusion is strict, even when one restricts our attention to differences 
of characteristic functions. For instance this test shows that one cannot go 
from the position of figure [TT] to only the central peg while the rational and 
integer linear tests are passed. This example is interesting in showing the 
impact of the board, for it is feasible in legal moves if we add to the english 
board the grey square on the upper right side. 



Figure 11: Impossible 

We present a smaller counterexample in figures [I2] and [I3] that enables 
easier direct computations. 



Figure 12: Starting position 



Figure 13: Ending position 



When g belongs to this intersection (i.e. the right-hand side of (fTTl) ) the 
denominators in a non-negative writing do not seem to be any worse than 
1/2. Here is the conjecture we make: 

Conjecture 6.1 If & has no is olated points, then 

y(6,Z)f|\/+(6,Q)ciV^+(6,Z). 
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Here is another related conjecture that may be easier to handle (and maybe 
easier to disprove!). 

Conjecture 6.2 Let B C be a basis ofV{&,Q). If & has no isolated 

points, then 

The condition on 6 cannot be removed since it is equivalent to V{&, Z) being 
of full rank. 

At this point, we have described the situation and we hope the reader is 
now able to understand properly what is what. The theory so far has two 
drawbacks: it draws only on properties of 1/ — Ij, and it does not use the 
order in which the moves are played. Our next criteria, the simple quadratic 
test, will not go beyond this abelian nature, but will break the first hurdle. 
It is better to investigate the game a bit further before exposing it. 

7 How integer linear programming is used 

The cone V^{&,Z) is determined by the set ^(6) of generators. Let us 
introduce the notation f for the function over ^((3) that is 1 in f and O 
everywhere else. We consider the map 

=F(^(6),Z) 

The integer linear program we write is simply to minimize any linear form 
of the {x{f))f subject to the constraints 

Vf e ^(e), x(f) > O, and *(F) = 1/ - Ij. 

The linear form we choose is usually '^fx{f) since we know what should be 
its value if a solution exists. 

8 Thickness of a move 

Given a problem, say from / to J, we define the thickness of the move f to 
be the maximum number of times this move can be used, whatever sequence 
of legal moves fi,f2, ■ ■ ■ ,fk we choose. This thickness is zero allover if the 



V(6,Z) 

E -(fif- 

fe^(6) 
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problem is not feasible. In general, given h G V^{&,J.), we shall speak of 
the thickness of f at h. Computing this quantity is naturally difficult, but 
we can bound it from above and even provide a uniform bound for it. The 
main Theorem reads as follows 

Theorem 8.1 Let h G V~^{&,Z), fo G ^(©) and n be a resource count on 
& such that (tt, fo) = 1. The move fo can appear at most (vr, h) in any writing 
of h as a linear combination of elements of ^(©) with non-negative integer 
coefficients. 

The scalar product (tt, /t) is defined in (l9|). We can derive absolute bounds 
from this Theorem by using a variant of a resource count already used by 
Conway. First note that we are interested only in the case h = Ij — 1 — J 
which implies that \h{A)\ < 1 for all A e &. Now let p = (a/5 - l)/2 be a 
solution of + z = 1. To each point (a, 6) G Z^, we associate the weight 
7r(a, 6) = pl'^l+l^l. Next, we drop our board & on in such a way that 
the middle point of fo be the (O, 0) element. The reader will check that the 
restriction of vr to (5 is a resource count on & which we denote again by tt. 
We have (tt, fo) = 1, while 

I (tt, h) I < (tt, le) < 8p + 13 = 17.944 ■ • ■ . 

This short argument show that the thickness of any move on any board 
is bounded above by 17. This is most probably a way too large majorant 
(reaching a thickness of 4 is already extremely difficult, and it can be shown 
on using better resource counts that the maximal thickness on the english 
board is at most 5), but it is universal, i.e. independant of the board we 
choose. 

A similar argument is also the main ingredient of [6j (see Theorem 3.1 
therein, with most probably a wrong computation at the end. The 26 of 
this result is to be replaced by a 34 but this leaves the rest of the argument 
intact), and is the basis on which rely the low complexity results. 

Given a problem, we can refine this upper bound by selecting a more 
appropriate resource count. Furthermore, once a majorant is given, say m, 
we can check whether 1/ — Ij — mfo is feasible or not (this means, whether 
it passes whichever test we select). If not, we decrement m and repeat the 
process. 

9 A simple quadratic test 

Let us consider the two following problems: we are to go from the left hand 
side position with only the black pegs (or with the grey peg added) to the 
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right hand side one with a sole black peg (or with the grey peg added). Both 
problems pass the positive integer test. The reader will easily check that the 




Figure 14: Starting position Figure 15: Ending position 

larger problem (with the grey peg) is in faet doable in legal moves, which 
implies that no test relying only on 1/ — Ij would be able to show the first 
problem to be impossible. The quadratic test we propose now is however 
able to show this impossibility. 

Let us start our description of the quadratic test. 



9.1 The geometricai support 

To each couple {A, B) & & x &, we associate a symbol A\S B, to which we 
add the property 

A\EiB ^ B\EiA. (13) 

We set 

e\se = {AmB,A,B eG}. (i4) 

We next consider functions on S Kl ©. We denote by A\S B the function 
that is 1 on A Kl 5 and O everywhere else. Note that A\ÉiB ^ B^A. We 
go from ^(6, Q)^ to ^(6 Kl 6, Q) by 

Kl : =^(e,Q) X ^(6,Q) ^ =^(6K16,Q) 

{91,92) ^ 91^92= Yl 9M)g2iB)A^ÉlB. 

A,Bee 

Notice that the value of 5-1 IS 5(2 on ^ Kl B is gi{A)g2{B) + gi{B)g2{A) if 
A 7^ 5 and gi{A)g2iA) ifA = B. 
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9.2 The efFect of legal moves 

Assume now that we can go from J to J by the legal move f G ^(©). We 
have 

1/ Kl 1, = (Ij + f) H (Ij + f) = Ij H Ij + f Kl Ij + Ij H f + f H f . 
On using the identity f Kl Ij = Ij Kl f , we reach 

1/ Kl 1/ = Ij Kl Ij + (21j + f) Kl f . 

We note that 



21j + f= J2 2^ + 



A&J 
KA)=0 



from which we infer 

1/ Kl 1/ = Ij Kl Ij + Ifl Kl f + 2i Kl f. (15) 

AeJ 
f(A)=o 

This is the equation we want to exploit; we do so in pretty much the same 
way we exploited ([T]). We set 

^(6 Kl 6) = {2i Kl f, A e ©,f e ^(6)/f(A) = 0} 

Udflet, fe ^(6)}. (16) 

Note that if f = P + Q - R then 

|f|Klf = PKlP + 2PKig + QKlQ-i?Kli?. (17) 
We define our cone by 

y+(6K16,Z)= J2 (18) 

cGfi?(6IX16) 

A problem being given by an initial position / and a final one J, the simple 
quadratic test consists is saying that 1/ Kl Ij — Ij Kl Ij G V'^{& Kl ©,Z), 
which can again be solved with integer linear programming. However the 
spaces are much larger, and the resolution becomes more troublesome. Note 
the following Lemma: 

Lemma 9.1 

|6 Kl 6| = |6|(|6| + l)/2 , |^(©K16)| = (|6| -2)|^(6)|. 
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Indeed, there are |i^((3)| moves of type |f| Kl f, and, for each f G ^(6), there 
are |6| — 3 moves of type Kl f with f{A) = 0. For the english board, the 
cardinality of |^((5 K &)\ is thus 2 356 for a board of 561 squares. 

We have already given an example showing that this test is sometimes 
better than the linear test with non-negative integer coefiicients but we show 
now that this is always the case. To do so, let us define 

{^o(<3Kie,Z) = ^Z-AKlA+ Z-2AK15 
w{eme,z) = ^ z-2A\sb. 

Then we can easily identify ^o{& K 6, Z)/W{& K 6, Z) with the space of 
integer valued functions on {A\Z1 A, A E &}, which we can in turn identify 
with &. By these identifications, we start with a function h G ^((3,Z), 
build h\Zi h E ^o{& Kl &, Z) and is next send to h. In particular, we get 

1/ K 1/ - Ij K Ij G 1/+(© K 6, Z) =^ 1/ - Ij G V^{&, Z). (19) 

The faet that this test is in faet strictly superior on some boards in shown 
by the problem described by figures [TH and [151 



10 A quadratic test, with flatness constraints 

If the simple quadratic test is stronger than the linear one with positive 
integers, it turns out when used to be lacking in efficiency. The last term 
in (ITSll can be written as 21k K f where K C & avoids the support of f. 
This is much better than saying that it is a linear combination of 2A K f, 
but it leads to 2l®l~^|^(©)| + |^(6)| generators! This is of course way too 
much and makes this new set of generators impractical. However, if JF is a 
succession of legal moves from I to J, we can write 

ljmij-lj\Silj = J2 <fM ^ f + E E ^f(^)2^ ^ ^ (20) 

f f A 

And we readily see that on this writing that the following inequalities are 
satisfied 

0<y^{A)<x{f). (21) 

We call them the flatness constraints. Despite their number, these constraints 
renders the quadratic test much more efficient. In faet, The z(f) are related 
to the usual linear moves by 

l,-lj = 5^x(f)f (22) 
f 
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(see the process that enabled us to prove (fTOl) ) and as such can be controlled 
in size by the thickness of f at 1/ — Ij, as defined in section[8l 

On an english board, the x(f)'s are seldomly larger than 4, and on arbi- 
trary board they are anyway bounded. 

Notice that if 1/ Kl 1/ — Ij Kl Ij passes this test, then actually, it can 
be written as a linear combination with non-negative integer coefficients of 
2/1 K f with f{A) = O and diagonal moves |f| K f. To realize such a writing, 
given f, simply collect together all A's for which yj{A) has a given value into 
a set A. Note that these sets A are not the same as the sets K we used at 
the very beginning of this section, but are of same use. 

The problem described by figures [16] and [17] goes through the quadratic 
test with no flatness constraints, but is shown impossible as soon as we add 
these constraints : 




Figure 16: Starting position Figure 17: Ending position 

This new test is the main novelty of this paper and is extremely effi- 
cient in practice, though it requires a processor to carry out the required 
computations. 

We end this part with three further examples of problems shown to be 
impossible via the quadratic test with flatness constraints. Here are two 
problems, with a same starting position but different ending positions. None 
of them go through the quadratic test with flatness constraints: 
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Figure 18: Starting 
position 



Figure 19: First 
ending position 



Figure 20: Second 
ending position 



The third example is to go from the initial position to the intermediate 
ending position. This is shown to be impossible via the quadratic test with 
flatness constraints, though it again passes the simple quadratic test. More- 
over, the problem to go from the initial position to the final ending position 
is feasible in legal moves. 



Figure 21: Starting 
position 



Figure 22: Interme- 
diate ending posi- 
tion 



Figure 23: Final 
ending position 



11 Additional constraints, a first draft 

Now that we have seen that the quadratic test with flatness constraints is so 
very efficient, it is tempting to try to add some further constraints. This is 
the topic of these two last sections, but this part is still very much in progress. 
The reader may get the impression that it is not so much in progress than 
more bluntly unfinished. After some months of efforts, I have not been able 
to derive a unifying setup for what look like protrusions of a hidden structure, 
which is why I deliver them in tnhis state. 
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The idea we follow is to add geometricai information to control as much 
as possible these new variables yf{A) in (l20l) . 
Let us start with a fundamental inequality. 

Proposition 11.1 Assume we can go from I to J in legal moves. Then 
there exists a writing af 1/ Kl 1/ — 1 j Kl 1 j ( as in (l20l ) ) such that for every 
A G & we have 

o<Y,yfi^)+ E ^(f) < 1^1 - l-^l- (23) 

f f/f(A)7^0 

See (l28l l and (l29ll for refinements. Let JF be a succession of legal moves from 
I to J. We set 

p{A,J^) = J2miA)+ Yl (24) 
where the ?/f(A)'s and the x(f)'s come from (l20ll . 

Proof: Given a move f , let us look at the situation of the board before using 
this move. There are four possibilities for A: 

• f{A) = 1, which means that A is on the board and participates to the 
move. It is counted in x(f) and nowhere else. 

• A is not on the board but is created by the move. It is counted in x{f) 
and nowhere else. 

• A is on the board but does not participate to the move. It is counted 
in y^{A) and nowhere else. 

• A is not on the board and not created by the move. It is not counted 
any where. 

The proof follows by using this remark and an induction on |/| — \J\. We 
have equality if and only if the last case above never occurs, which means 
that A is never absent from the position for two consecutive moves. o o o 

We have seen that we can have equality in (l23ll . but we can even show 
that the right hand side is on average of the correct order of magnitude. 
Indeed we have 

J^PiAJ^) = |/|-2 + |/|-3 + --- + |J|-l + 3(|/|-|J|) 

- (1,1 - ii|)^±M±l 
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since there are |/| — 2 points on the first move that are on the board but do 
not participate to the move, then |/| — 3, and so on. As a consequence 

I I y^gg I I 

This shows that (l23l) prevents too wide deviations from the mean, at least 
if |/| + |J| and |©| are of comparable size. We propose to improve on this 
double inequality in three ways. 



11.1 Using the speed at which a peg gets inside J 

We define the depth of the point A with respect to the position S containing 
it to be the minimum number Depth(yl, 5") of legal moves required to remove 
the peg in A. If A is not in S, we set Depth(74, S) = 0. Let us recall a 
classical Lemma. 

Lemma 11.1 (Leibniz) Ifthe sequence of legal moves fi, ^2, ■ ■ ■ ,fk goes from 
I to J , then the sequence of legal moves f^, . . . , f2, fi goes from & \J to &\I . 

It is enough to verify this property when k = 1 where it is obvious. Leibniz 
expressed this idea in a different manner: he started from the final position 
J and tried to recover the initial one by playing in reverse; he discovered it 
was the same game, provided one considered the empty squares as having 
a peg, and the ones with a peg as being empty. This is exactly what we 
shall consider. Indeed, given a point A out of our final position J, there is a 
minimal number a moves that will "bring" its peg inside J, or kill it, namely 
Depth(A,6 \ J). 

Let us select a minimal path from J to A. Its last move puts a peg in A, 
i.e. has A as point R since we could otherwise shorten this path. Moreover 
it does not use A anymore as point P or Q since we could again shorten the 
path. Consequently, for any A ^ J 

p(A,J^) < max(0, |/| - |J| -Depth(A(5 \ J) + !)• (25) 

If A is in J, we have Depth(y4, © \ J) = O so that (l23l) is stronger. 

Proof: Indeed A not in J implies Depth(A, 6 \ J) > 1. The Depth(A, 6 \ 
J) — 1 last moves cannot use A in any part of a move, hence we can use (l23ll 
with I J| + Depth(yl, G \ J) — 1 points as a final position instead of J if A 
is at some point of time on the board. Else, it is never here and the upper 
bound O is fine. o o o 
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We do not know of any precise mean of computing this depth, but we 
provide now a fast way to get an excellent lower bound. Let us consider the 
oriented graph & built on the set & and where we put an edge from A to B 
if there exists f e ^(©) such that f{A) = 1 and f{B) = —1. A minimal path 
that realizes Deptli(y4, &\ J) is readily transformed in a path from A to J 
on (5. Reciproquely from such a path from A to J on this graph, we deduce 
a position K by adding the required points P and Q necessary for the f s. 
The only problem is that this process may require to put several pegs on a 
same square (we do not have any example of such a situation). Denoting by 
S&{A, J) the distance on this graph, we have established that 

5&{A, J) < Depth(6 \ J) (26) 

Note that a final position L in case of 5© is reduced to a single point. The 
distance Sis{A, J) is now readily computed, by using the Dijkstra's algorithm 
for instance. 

Practically, to find a minorant of this depth, we proceed in two steps 
(with ^ = 6\ J): 

• We try every succession of 5 legal moves from S. 

• Concerning the remaining ones, we first build the set of points with 
Depth(A, S") < 5. If A G 5*5, we find the minimum of 5^{A,B) + 
Depth(i?, S) for every B E S^; this a first lower bound for Depth(A, S), 
but sometimes the lower bound 6 is simply better. 

11.2 Using the speed at which a point is reached by / 

Let us now examine the somewhat reciproqual situation, and try to get the 
minimum of legal moves from the set / that puts a peg in A. We need two 
pegs to create one, which means that the distance 6e,{A,I) is not a good 
lower bound anymore. We define the height Height(y4, J) of A with respect 
to / to be his minimal number, and set Height (A, I) = oo ii A can never be 
reached. Computing this Height(yl, J) is very difficult. 

Lemma 11.2 Let I be a subset of & and A be such that Height(y4, J) < oo. 
For any non-negative resource count n, we have (l/,7r) > tt^A). 

Proof: Indeed there is a set J which contains A and that is reachable from 
I. We thus have (1/, n) > (Ij, vr) which in turn is non less than vr(P) by the 
non-negativity assumption on vr. o o o 

Using Lemma 111.21 and some direct computations, we get the following 
height-diagram for the left-hand side position. 
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Figure 24: Starting po- 
sition 
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Figure 25: Height / - 
Depth 



We next provide an example on which Lemma ril.2l is not strong enough to 
decide whether some points have finite heights or not. This problem passes 
the linear integer test. We provide the height of each square (we simply 
computed all position attainable in 5 moves !). The two squares on the left- 
hand side (and the symmetric ones on the right-hand side) are rather clearly 
not reachable, but the test deduced from Lemma 111.21 fails to prove that. 
Even worse, we found for each of this square a position got from the first one 
in 5 moves and for which this square is not shown to be unreachable by this 
test. 
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Figure 26: Starting 
position 



Figure 27: Endin^ 
position 



Figure 28: Height 
-Depth 



Practically, to find a minorant of this height, we proceed in three steps: 

• We try every succession of 5 legal moves. 

• We use the Lemma 111.21 to determine those points that are guaran- 
teed to have infinite height. (We apply this test to all of the derived 
positions). 
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• Concerning the remaining ones, we first build the set of points with 
Height(v4, /) < 5. U A e & \ I5, we find the minimum of 6&{A, B) + 
Height(i?, /) for every B G /s; this a first lower bound for Height(A, J), 
but sometimes the lower bound 6 is simply better. 

Set 

"^{AJ.J) = mm(^\I\ - |J|,max(Depth(A© \ J) - 1,0) 

+ max(Height(A/) - 1,0)). (27) 

We have 

p{A,T) <\I\-\J\-'é'{A,I,J). (28) 

11.3 Using the speed at which a peg comes out of / 

We finally improve on the lower bound in (l23l) . The faet is that some points 
are so much within the starting position / that the peg on them cannot 
be eliminated before so many moves, and this is precisely how we defined 
Depth(A,/). We have then 

I)epth{A,I) <p{A,J^). (29) 

11.4 Final discussion 

We end this section with two remarks. First, both notions of depth and 
height use only one of the two positions of the problem, and this is a loss. 
For instance concerning height, if we manage to put a peg in a very far away 
square that is also far from our final position, it is probable that we shall not 
be able to bring it back to it; for instance, if the starting position is given by 
figure EHl it is likely that we cannot put a point in the lower left corner and 
finish as in figure [271 Secondly, constraints ( l28l ) and ( l29l l only avoid extremal 
cases, as we noted earlier, and there are only 2|(5| of them for a problem with 
about |(5p variables; in faet, if & has no isolated point. Lemma HH] yields 

\e\{\&\ - 2) < \^{& Kl 6)1 < 4|e|(|6| -2). 

This explains why these constraints are somewhat weak. 

12 Additional constraints 

Having in mind the counting argument displayed at the end of last section, 
we see that finding conditions on couples {A, A') of points would not increase 
too much the size of the problem but may yield more stringent constraints. 
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As of now, we have only found one such type of constraint, which applies 
to initial positions / such that 6 \ / is large enough. 

Let us start with some general considerations. Let Height(>l, A', I) be the 
minimum number of legal moves necessary to put a peg in each of A and A', 
starting from a board with pegs on all the points of /. We assign it value oo if 
no such succession exists. Note that the height-function does not behave like 
a distance, since we can have Height(74, A', I) > Height(A, 7) + Height(74', I). 
We formulate a conjecture: 

Conjecture 12.1 Reight{A, A' , I) > Height(A, /) + Height(Æ, /). 
A proof or disproof of this conjecture has sofar escaped the author. 

Lemma 12.1 Consider two points A and A' such that Height(A, A', /) = oo. 
Then 

p(A-^) + P(^',-^)- ^(f) < 1^1 - (30) 

f/f{A)f(A')^0 

Proof: Given a move f , let us look at the situation of the board before using 
this move. There are several cases: 

• A is on the board and is not moved by f. Then A' is not on the board, 
and may not be created by f. This move is counted in yf{A). 

• A' is on the board and is not moved by f. Then A is not on the board, 
and may not be created by f. This move is counted in yf{A'). 

• A is on the board and is moved by f. Then A' is not on the board and 
may be created. This move is counted in x{f). 

• A' is on the board and is moved by f. Then A is not on the board and 
may be created. This move is counted in x{j). 

ooo 

The question arises as to whether this Lemma leads or not to improvements, 
and we provide an example below showing that it indeed does. The geo- 
metricai faet that we have used is that a square can either contain a peg, 
or be empty, a fairly trivial information that was until now absent from our 
discussion. 

Before exposing our example, let us address rapidly the problem of com- 
puting couples {A,A') with }ieight{A, A' , I) — oo. 

Lemma 12.2 Let I be a subset of & and A and A' be two points of &. If 
there exists a non-negative resource countn, such that (l/jVr) < n{A)+7:{A'), 
then Height(/1, A', I) = oo. 



25 



We can improve on this criteria: simply form all positions derived from / 
by one (or any fixed number) legal move, and apply this criteria to each of 
them. 

Here is a problem that is shown impossible by using this criteria, though 
it passes the quadratic integer test with flatness constraints: 



1 -1 
-1 1 



-1 -1 
1 1 



Figure 29: Starting 
position 



Figure 30: Ending 
position 



Figure 31: Height 
-Depth 



This example is also interesting because of the square with an interroga- 
tion dot: it is "clearly" of infinite height, but our automatic process is not 
able to conclude. Here is the list of couples with Height (A, A', I) = oo that 
we have found: 



A 


A' 


(3,1) 


(1,4) 


, (2, 3), (3, 2), (3, 3) 


(4,1) 


(1,4) 


\ {1, 5), {2, 3), {2, 4), (2, 5), (3, 1), (3, 3), (3, 4), (3, 6), 




(3,7) 


,(4, 2), (4, 3), (4, 4), (4, 5), (4, 6), (4, 7) 


(5,1) 


(1,4) 


, (2, 3), (3, 1), (3, 2), (3, 3), (4, 1), (4, 2), (5, 2), (5, 3) 


(5,2) 


(3,1) 


, (4, 1) 


(5,3) 


(3,1) 




(5,4) 


(4,1) 




(6,3) 


(3,1) 


,(4,1), (5,1) 


(6,4) 


(4,1) 


,(5,1) 


(6,5) 


(4,1) 




(7,3) 


(1,4) 


, (2, 3), (3, 1), (3, 2), (4, 1), (4, 2), (4, 3), (4, 6), (5, 1), 




(5,2) 


,(5, 3), (6, 3), (6, 4), (6, 5), (7, 4), (7, 5) 


(7,4) 


(1,4) 


, (2, 4), (3, 1), (4, 1), (4, 2), (4, 4), (4, 6), (4, 7), (5, 1), 




(5,2) 


,(5, 4), (5, 6), (6, 3), (6, 4), (7, 5) 


(7,5) 


(4,1) 


,(5,1), (6, 3), (6, 5) 



Lemma [l2.1l is of course of fairly limited use: we need the starting position 
to leave free enough squares on the board. However, it shows how more 
geometricai arguments may be used to get improvements! Our journey ends 
here. 
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